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Abstract. We introduce a general class of models for charaterizing the non-Gaussian 
properties of foreground contaminants in the cosmic microwave background with view towards 
the removal of the non-primordial non-Gaussian signal from the primordial one. This is 
important not only for treating temperature maps but also for characterizing the nature and 
origin of the primordial cosmological perturbations and thus establishing a theory of the early 
universe. 



1. Introduction 

Cosmology is presently a very active field because of the large number of observations that are 
becoming available and that will allow us to characterize with great precision the nature and 
physical origin of the primordial cosmological perturbations. A key question is whether these 
primordial perturbations were Gaussian. 

The currently favoured best-fit cosmological models, as supported by the recent Wilkinson 
Microwave Anisotropy Probe (WMAP) data, are in agreement with the primordial cosmological 
perturbations being Gaussian to great accuracy pQ. With the WMAP data there has been a 
claimed detection of primordial non-Gaussianity that is currently disputed [2]. However, even 
the conservative official analysis hints at something, the value of /jvx lying within 1.9cr from the 
null value. 

Either a detection or a more stringent constraint of /jvl would have profound implications for 
our understanding of the physics of the early Universe. In particular, it would be an observational 
discriminant among competing models for the generation of primordial fluctuations. An /tvl ^ 
would be a very interesting challenge for cosmology since it would require an extensive revision 
of the standard cosmological model. Present observations have ruled out a great number of 
cosmological models already. All currently favoured models would be ruled out by such a 
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detection. A non-detection of //vx with an improved constraint on the allowed values would 
strongly restrict the fluctuation generation mechanisms. 

Due to the great improvement in sensitivity of the European Space Agency's Planck satellite 
over WMAP, understanding the nature and in particular non-Gaussian properties of foregrounds 
will be indispensable for the Planck /jvl analysis. For WMAP, foregrounds were not an issue 
because no signal was detected. Currently estimators exist [3, 4j that are proven optimal in 
the absence of foreground contamination and much work has been devoted to their efficient 
implementation [S]. However, to date there has been no comprehensive study of the impact 
of foregrounds. Foreground sources distributed in an anisotropic way are susceptible of being 
misidentified as primordial non-Gaussianity. In order to extract from the data an accurate and 
significant measure of Jnl, new estimators are needed that divide the data and are designed 
to mask possible contaminants. First, however, it is required to characterize the statistical 
properties of these spurious sources of non-Gaussianity, which requires understanding the non- 
linear properties of foregrounds. 

This contribution reports on work in progress where we formulate a general class of models for 
characterizing the non-Gaussian properties of foreground contaminants in the cosmic microwave 
background (CMB) [BJ. 

2. Measuring non-Gaussianity 

In linear perturbation theory, the temperature anisotropy of the CMB is derived from metric 
perturbations by the relation AT/T oc gxA&, where the radiation transfer function gx depends 
on the scale considered!! Linear perturbation theory is a valid approximation because the physics 
of the pertubations generated during inflation is extremely weakly coupled. 

In the simplest models of inflation, metric perturbations are produced by quantum 
fluctuations of the inflaton field. If the physics that governs the evolution of the metric 
perturbations is linear, since quantum fluctuations are Gaussian then so will temperature 
fluctuations be Gaussian. However, non-linear coupling between long and short wavelength 
fluctuations of the inflaton can generate weak non-Gaussianities. These will propagate into the 
metric perturbations and consequently into the temperature anisotropy. The primordial metric 
perturbations $> can be characterized approximately by a non-linear coupling parameter, /atl, 
such that [7] 



where $l(x) are Gaussian linear perturbations with ($x(x)) = 0. In many models (e.g. single- 
field slow roll inflation [HJ [U]), the non- linear corrections are too small to be detected. Other 
models (e.g. multiple field inflation [101 lllj ) can generate stronger and potentially detectable 
non-linearity. A Gaussian distribution is characterized by vanishing odd-order moments and 
even-order moments defined only by the variance. Non-Gaussianity can thus be measured 
by correlations among an odd number of points. Here we shall use the angular three-point 
correlation function or bispectrum as our statistical tool since it is sensitive to weakly non- 
Gaussian fluctuations. 

Expanding the temperature anisotropy in momentum space 



For temperature fuctuations on super-horizon scales the Sachs-Wolfe efect dominates and <?t is a function of 
the parameter ui of the equation of state. At the decoupling epoch it reduces to gr = —1/3 in the case of adiabatic 
fluctuations. On sub-horizon scales, however, gr oscillates and we need to solve the Boltzmann photon transport 
equations coupled to the Einstein equations for gr, which is found to have as parameter the decoupling time. 
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where the harmonic coefficients are 

aim = M^Y J ^$(k) 5T; (k)F4(0), (3) 
the bispectrum is defined as 
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Here 

G ™™ = / dnY hmi (n)Y l2m2 (Q)Y hm3 (n) (5) 

is the Gaunt integral and 

6^/2/3 = (^j J dx dki dk 2 dk 3 (xkik 2 k 3 ) 2 B< s> (ki,k 2 ,k 3 ) 
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is the reduced bispectrum defined in terms of the radiation transfer functions gT n the Bessel 
functions of fractional order ji and the primordial bispectrum B$(ki,k 2 ,k 3 ). In the flat-sky 
approximation, which amounts to (l,m) — > £, where £ is a two-dimensional wave vector in the 
sky, the bispectrum reduces to 

(a(£i)a(£ 2 )a(l 3 )) = (2tt) 2 5 2 (£ 1 + £ 2 + £ 3 )B(£ U £ 2 ,£ 3 ) (7) 

which is non-vanishing only for closed triangle configurations in momentum space. Since in this 
approximation G™^™ 3 -» (2n) 2 5 2 (£ 1 + £ 2 + £ 3 ), then -» B(£ u £ 2 ,£ 3 ) p]. 

The signal-to-noise squared for the bispetrum gives us an indication of the dominant 
triangle configuration to the non-Gaussian signal. We assume a scale invariant 
power spectrum P(£) oc l~ 2 in the flat sky approximation and B{£\,l 2: i 3 ) = 
Inl [P{h)P{li) +P{£2)P{h) +P(4)P(4)] as a result of the local non-Gaussianity ansatz. 
Then 
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For squeezed triangles, where the length of the two large sides is of order t ma x and that of the 
smaller side £ m i n , we find that 
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For equilateral triangles with the length of the sides of order 
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Comparing the two results we observe that the dominant contribution to the signal comes from 
the squeezed triangles, which introduce the logarithmic boost factor and thus couple the small- 
scale with the large-scale anisotropics. Should only equilateral triangles be considered, then the 
description of the signal would be equivalent to the separate description from different patches 
in the sky, which would mean losing most of the signal. The squeezed triangle configurations 
correspond to the local case where the non-Gaussianity was created primarily on super-horizon 
scales, as opposed to the equilateral case where the non-Gaussianity is primarily created at 
horizon-crossing. Qualitatively an /jvx 7^ can be understood as a modulation of the small- 
scale power (on scales near the resolution of the map) by the large-scale signal. Exploring further 
this idea, we can construct a cubic estimator of f^L by filtering the signal into two contributions, 
namely T\ ow £ and Thigh e> an d considering the integral 

E = j dVL T low £Th ig h £• (11) 

This is but a caricature of the optimal estimator given in Ref. [13] with the ideal weighting, which 
here serves illustration purposes. The large-scale variation in T^ igh e reflects spatial dependence 
in the power spectrum P(£, O,) where I is large and the variation in f2 is slow. In a Gaussian 
model, such perceived spatial dependence is simply noise from the cosmic variance and thus 
fortuituous. In a model with /jvl ^ 0, however, it is correlated with the large-scale temperature 
anisotropy. 

Astronomical microwave sources as well as instrumental effects can produce spurious non- 
Gaussian signals which would be measured together with the primordial signal by the bispetrum 
and thus susceptible of being confused for primordial non-Gaussianity. Examples of such sources 
are [Hj: the synchroton emmision, which results from the acceleration of cosmic ray electrons in 
magnetic fields; the free-free emission, which is produced by electron-ion scattering and is thus 
correlated with Hydrogen emission lines; and thermal galactic dust emission. 

The separation of the CMB and foreground signal components has relied on their differing 
spectral and spacial distribution. However, correlations at different scales of the power of 
non-primordial foreground signals could not be separated from the CMB by the canonical 
methods. Here we aim to show that a detectable non-Gaussian signal could be mimicked by 
such correlations. 

3. Faking non-Gaussianity 

Since the a; m 's are random quantities for the same sky realization, we must devise a way of 
generating maps of the sky with built-in non-linear foregrounds. We use motifs of a fixed shape 
and distributed in a probabilistic manner in the sky to model foreground sources. This resembles 
the galaxy counting method for the study of the large-scale structure formation |15j . 

Foreground realizations are generated in two steps. First a point £ for the centre of a motif 
is randomly selected according to a probability density function Q(£) in the simplest case in 
an uncorrelated Poissonian manner. Then to each point 9 we assign a conditional probability 
P(6,£) of finding the point for a motif centered at This probability is independent of how 
many points are already inside the motif. We assign the value one to the motif interiors and 
zero to the regions between motifs. Assuming that the foregrounds are optically thin, i.e. that 
they do not mask the effects of each other, we take overlapping motifs to add. We also need 
to properly normalize these quantities so as to relate to the foreground maps. Measuring the 
mean number density of motifs, the probability density Q(£) is such that when integrated over 
the entire sky it yields the average number of motifs (N) 



d 2 £, Q{£) = (number of motifs) 



(12) 



The conditional probability P(6,£) yields the temperature at a point given that a motif 
is centred at Thus correlations among the £'s as encoded in Q will model the large-scale 
structure, while correlations among the #'s in the same motif as encoded in P will model the 
small-scale structure. 

We calculate the two and three-point correlation functions of the statistical map in a 
perturbative manner using Feynman diagrams. These correlation functions are interpreted 
respectively as the power spectrum and the bispetrum by assigning a temperature to the motifs. 
This temperature is contrained: from below by the minimum temperature that a foreground 
distribution must have to account for the remaining signal after removal of the astrophysical 
sources described previously; and from above by the maximum value of the power spectrum that 
a spurious source could have so as not to overshadow the background radiation. 

We define the two-point correlation function of the final map V(0i,0 2 ) by the sum over the 
two contributions to the probability that Q\ and 02 are within a motif 



where 



V{Oi, 6 2 ) = pi,i(6>i, 6> 2 ) + p 2 (0i, 2 ), (13) 

pi,i = d 2 ^J d% P(e 1 ,^ 1 )Q^ 1 )p(e 2 ^ 2 )Q(^) 

p2 = h. j ^i p (^i) p ( 2,£i)Q(£i). (14) 

Here p^i (0i, 02) is the probability that the points 0i and 6 2 are within different motifs centred at 
£i and £ 2 , an d P2(0i, #2) is the probility that the points are within the same motif. Analogously 
we define the three-point correlation function B(0\, 02, 0%) by the sum of the three contributions 

B(01, 2 , 3 ) = Pl,l,l(01, 02, 3 ) + P2,l(01, 02, 03) + P 3 (01, 02, 03) (15) 
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are defined analogously. The combinatorial factors preceding the integrals are the number 
of equivalent motif configuration for indistinguishable patches in the £-space which must be 
introduced to avoid multiple counting. We introduce a factor 1/uq ! where uq is the number of 
Q's, and a factor 1/np! for each £ where np is the number of P's shared by the same £. For 
convenience we also define the one-point function 

Pi(0i) = J d% P(0i,^)Q(^). (17) 



4. A Simple Model: Uniform Distribution 

For concreteness we assume the motifs to be circles of radius 9 c i rc , with sharp edges P(0, £) = 
NpQ[9 c i rc — |£ — 0|] and uniformly distributed in the sky Q(£) = const = Nq. Thus Nq = 

(N) /A sky and A^ P 

— P-motifi which in this case are constant inside the motif and independent 
of the position of the motifs in the sky. [To be dimensionally correct, we must have T mot if 
normalized to the average temperature of the sky T s k y ] We will also discuss a non- uniform 



distribution of the motifs in the sky, motivated by the variation of the distribution of the large- 
scale structure along the latitude. This will be the working case for introducing dependence of 
the temperature on the position of the motifs as encoded in the distribution probability. [For a 
refinement of the functional form see Ref. [6].] 

We compute the two and three-point correlator functions in real space for the statistical 
ensemble just described. The one-point function becomes 

Pi(0l) = N P N Q J d 2 x x Q[9 circ - \xi\] = N P N Q A motif (18) 

where A mo uf = it9 2 irc . For the power spectrum we must calculate pi 5 i and p 2 , which can be 
written as 

Pl,l(0i,0 2 ) = \ { N 2 P N 2 Q j d 2 x x Q[9 arc - \ Xl \] J d 2 x 2 Q[e circ - \x 2 \] = ^p? (19) 

p 2 (6>i,6> 2 ) = ^NpN Q J d 2 x 1 e[e circ -\x 1 \]e[9 circ -\x 1 + e 12 \]. (20) 

Here Xi = £j — 6{ and Gij = Oi — Oj. For the bispectrum we observe that both pi,i,i and p 2 ,i 
can be expressed in terms of pi and p 2 as follows 

Pi,i,i(0i,02,0 3 ) = ^NpN^ j d 2 Xl 9[^ rc - |aji|] J d 2 x 2 e[e arc - \x 2 \] 
x J d 2 x 3 Q[9 circ - \x 3 \] 
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(21) 



p2,i(0i,e 2 ,B 3 ) = ^ n p n q J d2x i ®[0drc-\xi\] &[e ci rc- 1*1 +e 12 \] 

x J d 2 x 3 Q[9 c i rc — \x 3 \] + (two permutations) 
Pi [p 2 (6>i, 6 2 ) + p 2 (0 2 , 6> 3 ) + p 2 (6> 3 , Oi)] , (22) 
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so we have in addition to compute p3 only 



P3 



(Gi,e 2 ,0 3 ) = ^NpN Q J d 2 x 1 @[9 circ -\x 1 \] Q[9 circ - \x x + 12 \] Q[9 circ - \x x + 6» 13 |]. (23) 



These quantities describe the probability that one, two or three points, denoted by G x , 62 and 
3 , are within the same motif centred at points £'s which are distributed in the sky according to 
Q. These probabilities can also be interpreted as the measure of the overlap of motifs centred at 
each of these points #'s. The integrals of the step functions about the distance from the centres 
at £'s can thus be visualized as the area of the intersection of motifs centred at the #'s. [For 
more details see Ref. [6].] We find that 

Pi,l(0i,0 2 ) = ^P? (24) 
p 2 (6>i,0 2 ) = ^NpN Q A ln2 (25) 
where ^4in 2 is the area of intersection between the motifs centred at 9 X and # 2 

^4in 2 = I 20Ev„ arctan 
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These contribute to the power spectrum. Moreover, we find that 

Pl,l,l(01, 02,03) = T^P? (27) 
P2,l(01,02,0 3 ) = ^Pl[p2(01,0 2 )+P2(02,0 3 )+P2(03,0l)] (28) 

p 3 (0i,02,0 3 ) = 7^iVpiVQ^in2n3 (29) 

where ^4in2n3 is the area of intersection among the three motifs 
1 



^4in2n3 



, v (|012| 2 + |0 23 | 2 + |013| 2 ) 2 - 2 (|0l 2 | 4 + |0 23 | 4 + |013| 4 ) - 2 A motif 
X0 [26 c i rc — |012|] © [2^circ ~ | ^23 1] © [26circ ~ |^13 1] 
+ 2^ ln2 ® fflcirc — |023|] © [26circ ~ Wl3\\ 
+ 2^ 2n3 ® fflcirc — |012|] © [29circ ~ |^13 1] 

+ ^in30[2^ c -|0i 2 |]0[2^ rc -|023|]- (30) 
These contribute to the bispectrum. It follows that 

V(6 1 ,e 2 ) = lp 2 + p 2 (0 1 ,0 2 ) (31) 

#(01,02,03) = ^Pl + ^P1[P2(01,02)+P 2 (02,03)+P2(03,0l)]+P3(01,02,0 3 )- (32) 

We now proceed to compute the value of f^L produced. The parameter f^L characterizes 
the amplitude of the temperature non-Gaussianity since it couples to the quadratic term of 
the expansion of the temperature fluctuation T = AT/T about a Gaussian distribution Tl. 
Assuming the Sachs-Wolfe approximation (AT)/T = — (l/3)A<I?/c 2 on all scales and for an 
infinitely thin surface of last scattering, we have in real space that 

T(0) = T L (0) + 3f NL [rl{0) - (T L (0)) 2 ] . (33) 

Here the Gaussian distribution has zero mean, (Tl{6)) = 0, from which it follows to leading 
order that 

2 ) = (T(e 1 )T(e 2 )) = (r L («i)T L (e 2 )) . (34) 

Since a Gaussian distribution has vanishing odd-order momenta, we find for the three-point 
correlation function that 

B(0i,02,0 3 ) = <T(0i)T(0 2 )T(0 3 )> 
= GfNL [V(0i, 03)^(03, 2 ) + V(02, 0i)V(0i,e 3 ) + V(0 3 , 2 )V(02, 0i)] • (35) 

We distinguish three cases for the three possible relations among the distances between the points 
9's, namely 1) #12, #13, #23 > 2# c j rc ; 2) 9±2 < 26 c i rc and #13,6*23 > 26> c j rc ; 3) #12, #13, #23 < 26> c j rc . 
Let us fix the size of the motif and for concreteness take the radius to be of the size of the 
resolution of the map, i.e. 9 c i rc ~ 10~ 6 . 



4-1- Case #12, #13, #23 > 2Qcirc 

This is the case where each point 0j is within a different motif and consequently only the 
one-point function pi and the contributions to the two and three-point functions which can be 
expressed in terms of pi, i.e. pi 5 i and pi,i,i, do not vanish. We find that 

p(0i,e 2 ) = r(o 1 ,o 3 ) = r(o 2 ,o 3 ) = p 1 , 1 = ^pl 

B(0i,02,03) = Pi,i,i = ^P?, ( 36 ) 
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from which it follows that 



/JVL "6 3^(01, 2 )P(02, 03) "27pi' Pl>L (J7j 
^.2. Case #12 < 29 circ and #13, #23 > 2# circ 

This is the case where two points, here for concreteness 0i and 02, are within the same motif. In 
addition to the one-point function and the contributions derived from it, we have p2(0i, 02) 7^ 
and p2,i(0i, 02, 03) 7^ from the contribution to the three-point function of the two-point 
function between 6\ and 02 only. Hence 

v(e 1 ,e 2 ) = Pi,i+p 2 (0i,0 2 ) = ^p? + p2, r(e 1 ,e 3 ) = r(e 2 ,e 3 ) = p lil = ^pl 

B(0!,0 2 ,0 3 ) = Pl ! l,l+P 2 ,l(01,02,03) = ^P? + ^PlP 2 (01,02), (38) 

and 

f = 1 g(0i,02,0 3 ) = l ipf + ^pip 2 (0i,0 2 ) 

lNL 62V(e 1 ,e 2 )v(e 2 ,e 3 ) + r(e 1 ,e 3 )r(e 2 ,e 3 ) 6 2(±p? + p 2 (o 1 ,o 2 )) \p\ + \pf 1 ) 



We note that pf oc NqA moti f and p2 oc NqA\^ 2 < NQA moti f, where Nq is the mean number 
density of motifs in the sky. The relative magnitude of pf and p2 depends on the relation between 
the mean number of motifs (N) and the area of the motif A mot if. For convenience we define 
a = NQA mot if = (N) A mot if /A s ky. We must further distinguish between the following two cases. 
If a < a 2 , i.e. (N) > A s k y /A mot if ~ 10 13 , then p 2 < a < p\. This suggests that a minimum 
density of motifs is required in order for the one-point function to dominate, i.e. in order for two 
points which are close enough to be within the same motif to be also within a second, and thus 
necessarily overlapping, motif. On the other hand if a > a 2 , i.e. (N) < A s k y / A mot if ~ 10 13 , 
then both pf,p2 < Here in order to discriminate the relative magnitude, we need in addition 
to precise the relation between -Am2 and A moti f, which depends on #12. If #12 is sufficiently small 
so that ^4in2 ~ A moti f, then p2 ~ a and consequently p2 > p\- However, if #12 is very close to 
2# c iro then A\^ 2 <C A moti f and consequently p2 <C a, so that if p2 < a 2 then p2 < p 2 . This is 
equivalent to the case a < a 2 . 

Hence, if (N) > 10 13 , or (N) < 10 13 and # i2 ~ 2# circ , then 

'■"■~l$ = hk Pl>1 - (40) 



whereas if (TV) < 10 13 and #12 < 29 circ , then 

1 ^Pip 2 (0i,0 
6 p 2 (0i,0 2 )p 2 12 pi 



l ^PiP 2 (0i,0 2 ) 11 
JUL ~ T~~T7i TTyZT = To I - ' Pi < I" ( 4i ) 



4-3. Case 9i2,6l3,023 < ^dcirc 

This is the case where the three points are within the same motif. All the terms contribute to 
the two and three-point functions. In order to determine which contribution dominates for a 
mean number of motifs (N) distributed in the sky according to Q, we do an analysis similar to 
that above. Thus if (N) > 10 13 



fm ~\$r*h P1>L (42) 



However, if (N) < 10 13 we find for 12 < 6* 13 , 6» 23 that 

1 P 3(fli,fl 2 ,fl 3 ) = 1 ~\A motlS + \Mrq 
JNL ~ 6 2p 2 (6>i,0 2 )p 2 (0 2 ,6> 3 ) 18 N P N Q A ln2 A 2n3 ' { ' 

while for 0\ 2 ~ 623 ~ #13 we find that 

1 p 3 (0i,fl 2 ,fl3) 1 &6\ 2 - \A motl} + f A ln2 

lNL ~ 6 3p2(0i,e 2 )p 2 (e 2 ,e 3 ) 27 N P N Q Aj n2 ' 1 > 

For small 0\ 2 the term {\f^Ji)0\ 2 is subdominant, so for both cases we need only compare 
A\^ 2 with A mot if. We present the results for the former case, with those for the latter following 
straighfowardly. For < eq , where eq is the distance such that (3/2)^4i n2 = (1/2) A mot if , we 
find that (3/2)^4i n2 > (l/2)A mot if and hence that 

f NL ~ 1 ^ = 1 I . (45) 

JNL 18 N P N Q A ln2 A 2m 18 NpNqAm V ; 

For eg < ^i 2 < 20 c i rc the contribution of (^3/4)^ becomes more important but the dominant 
contribution is that of (l/2)A mot if . Moreover since ^4in 2 is a decreasing function of 0\ 2 except 
for 0\ 2 very close to 20 c i rc , then ^4in2^4 2 n3 > ^in3 anc ^ consequently — l/Ai n2 A 2n 3 > — l/^4 2n3 . 
Thus 

f > ^ ~\A-motif 1 —\ Amotif 

JNL ~ 18 N P N Q A ln2 A 2m 18 NpN Q A 2 2n3 { ) 

A similar calculation would follow for the case 0\ 2 ~ 2 3 ~ ^i 3 . 



Since these distances must be such that < #i 2 ,#i 3 ,# 23 < 20 c i rc , then — oo < Jnl < oo, i.e. 
we can generate arbitrary non-Gaussianity. We need, however, to be able to constrain /jvl to a 
finite interval. This requires the use of the foreground maps to constrain the parameters for the 
case of the uniform model, as well as the functional form of the probability functions for more 
realistic models. 



5. Constraining the Model 

Constraints on the parameters of the model can be extracted from the properties of foreground 
maps. Spatial templates from the WMAP data produced for synchroton, free-free and dust 
emission show a temperature distribution strongly dependent on the latitute. In order to account 
for this observation we must consider a non-uniform probability density that can capture the 
qualitative aspects of these dependence. The simplest case of a non-uniform distribution of 



motifs in the sky is that where Q{6) is a slowly varying function acroos the scale of the motif, 
so that the probabilities are changed to 



and to 

Pl,l,l(01, 02> ^3) 
P2,l(#l, 02, #3) 

p 3 (Oi,e 2 ,e 3 ) 

The changes in the two and three-point correlation functions, as well as in /nl, are 
straightfoward. Assuming that the contribution of the overlapping motifs add in the intersection 
region (what we called the "optically thin" approximation), then the value of Q in the mid 
distance between the centres of two overlapping motifs is the sum of the value of Q at the centre 
of each motif, i.e. 

=G(0i)+Q(0j)- (52) 
Since the scale on which Q (£) varies is large compared to A mot if, then 

m = ZM > (53) 

J- motif ^sky 

so that the same normalization condition holds. Possible functional forms for Q (£) will be 
explored in Ref. [BJ. For the purpose of this section, however, we will not need to be more 
specific. The conditional probability becomes P(0,£) = T(0)@ [9 c i rc — |£ — 0|] which will be 
simply the average temperature at the point £ where the motif on which lies is centred. 

By quantifying the clean regions of these maps, we can also set bounds on the foreground 
temperature. From a carefull examination of the five-year data maps [5] we extracted the 
minimum temperature To that a foreground source should have in order to saturate the map with 
as close to a single temperature value as possible. Taking the minimum value of the temperature 
of the three foregrounds for each of the five frequency bands analysed, we found by eye inspection 
that saturation of the maps was achieved for To a few \xK. Hence the temperature at each point 
of the motif distribution must be such that T(£) > To. On the other hand, the temperature of a 
single spot must be highly constrained so as not to stand out. This observation is captured by 
the condition that the power spectrum generated by the motif distribution does not outshine 
the background radiation. If the motif distribution produced a large power spectrum, as that 
produced by small and bright sources, then non-Gaussianity would be obvious. We thus set 

p(e it 0j) < (T /T sky ) 2 . 

The three cases discriminated in the previous section will now be combined with the 
constraints on the temperature obtained from the foreground maps. Apart from factors of order 
0(1), the results will be the same for both a uniform and a non-uniform probability distribution 
of motifs in the sky. For simplicity we will analyse the results from the uniform distribution. 
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5.1. Case #12, #13, #23 > 29 circ 

Here for any two pair of motif centres i,j, the constraints on the temperature yield that 

V(9 i} Oj) = \N 2 P N 2 Q tf motlf > ~ (J^j A 2 motlf N 2 Q (54) 



and simultaneously that 



P!0,.e,i<:(^) 2 (55) 

L sky J 



1 / A \ 2 

/ at\2 / Si-motif \ 



from which it follows that 

iWHT^) <L (56) 

V sky / 

This relation constrains the size of the motifs given their average number (N) distributed on 
the sky. Using the relations above in the expression for /jvl we find that 

111 1 

fNL ~ 27 N P A motif N Q < 27 (N) (A motif /A sky )(T /T sky ) ' (5 } 

Given the constraint from the power spectrum, we find that the magnitude of Jml is further 
determined by (To/T s k y ) only. Here /jvx is thus a measure of the deviation of To from the 
average sky temperature. 

5.2. Case #12 < 29 c i rc and 6*13,6*23 > 20drc 

Here for (JV> > 10 13 , or (N) < 10 13 and 12 ~ 29 circ 



nO^O^-^PO^O,) p i 9 2 .0 :ir =±N 2 Al otlf N 2 2>U-^ Al^jN 2 . CN) 

\ sky J 



which reduces to the previous case. For (N) < 10 13 and #12 «C Qdrc however, 
P(G 1 ,e 2 ) ~ \N 2 P A ln2 N Q > (^-) A ln2 N Q 

^ \J- sky J 



(59) 



which together with Eqn. ([55]) yields 



\(N) ^ < 1 (60) 

^ -"-sky 

This relation constrains the overlapping between any two motifs given the average number of 
motifs. Then for /jvl we find that 

f NT „ 1 < . (61) 

12 N P A moUf N Q 12 {N){A motif /A sky )(T /T sky ) v ; 

This relation sets a constraint on fw l of the same order of magnitude as the previous case. 



5.3. Case #12, #13 > #23 < 29drc 

Here for (N) > 10 13 we find the same result as in the first case where the constraint on the 
power spectrum is given by Eqn. (I56p . For (N) < 10 13 the constraint on the power spectrum 
is the same as in the second example of the second case and given by Eqn. (I60p . We further 
discriminate between two cases, namely #12 < #13, #23 and #12 ~ #13 ~ #23 > and each case for 
two regimes of the parameter 9\2- Thus for 612 < #13, #23 we find that 

, 1 ^4in2n3 J -\A- mot jf + \Mpg , , 

JNL ~ 18 N P A ln2 A 2n3 N Q 18 (N) A ln2 (A 2n3 /A sky )(T Q /T sky ) ' 1 j 



which for 612 < eq becomes 



fNL < 18 (N) (A 2n3 /A sky )(T /T sky ) ' (63) 



whereas for 9 eq < 612 < 20 c irc 

f > —2~A moti j 1 —^Amotif 

JNL * 18 (N) A ln2 (A 2n3 /A sky ){T /T sky ) 18 (N) A 2n3 (A 2n3 /A sky )(T /T sky ) ' 1 ; 

A similar calculation would follow straightfowardly for the case 612 ~ #13 ~ ^23- Eqn. (I63p sets 
a weaker constraint on /nl than Eqn. (jSTI) by the order of magnitude of ^2n3 /^matif- Since 
^4in2 is predominantly a decreasing function of #12, Eqn. (|64|) sets a weaker still constrain while 
allowing for a negative correlation among the three motifs. 



6. Discussion 

We propose a simple family of models for mimicking foreground sources that could contaminate 
the non-Gaussian signal of the CMB. This contamination could lead to the misidentification 
of a detection of a non-Gaussian signal for primordial when in fact we would be looking at 
the spurious signal from late-time, non-linear sources. Qualitatively non-gaussian aspects of 
foregrounds that are likely to give a significant signal of Jnl 7^ result from the modulation of 
the small-scale power by the large-scale power. 

Our model allows to generate foreground maps by distributing motifs in the sky according to 
a probability density Q(£) and correlating them according to a conditional probability P(9, of 
finding one, two or three points inside the same motif. The statistical properties of the resulting 
maps are determined by the two and three-point correlation functions, which are calculated 
by simply evaluated tree- level Feynman diagrams. We find the expression for /jvx in terms of 
the parameters of Q(£) and P(6,$,), namely the mean number of motifs and the intersection 
areas among two and three motifs. We suggest a prescription for introducing temperature in 
order to interpret the statistical properties of the motif ensemble as statistical properties of the 
temperature anisotropics. We also indicate how to use the foreground maps to constrain the 
normalization factors and accordingly we constrained the values for fNL- 

In the forthcoming paper we will use the model to generate concrete mock foregroun maps 
consistent with the level of foreground contamination observed and make detailed analysis of 
the impact on /jvx for Planck and other experiments. 
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